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Abstract 

We investigate the effects of spin fluctuations on the tunneling spectra of the NS 
junction. In the high junction resistance limit, the dip/hump structure observed in 
ARPES data for the high Tc superconductors is reproduced in a RPA treatment 
of the t-t'-J model. It is shown that the dip/hump structure weakens as doping 
increases as reflected in the data. In the other limit, we predict that the zero bias 
Andreev peak can coexist with the dip/hump structure. Furthermore, the c-axis 
tunneling spectra is found to be very similar to a recent STM data once these 
fluctuations are included. 

{PACS numbers: 74.20.-z, 74.50. +r, 74.80.FP, 74.20.Mn } 

While a great details about the spectral function for high-T c superconductors (HTS) have 
been revealed by the angle resolved photoemission spectroscopy (ARPES) Q, a complete de- 
scription of the superconducting state requires knowledge of the anomalous Green's function. 
Conventionally the tunneling spectroscopy has been considered as one of the tools which can 
probe the anomalous Green's function. In particular, measuring the subgap conductance of a 
junction consisting of a normal metal and a superconductor (NS) is the most convenient con- 
figuration for such a purpose. To calculate the conductance in the NS configuration, Blonder, 
Tinkham, and Klapwijk[^](BTK) developed a formalism using Bogoliubov-de Gennes (BdG) 
mean-field equations. The BTK theory has been phenomenologically extended to investigate 
the tunneling phenomena in various NS junctions^. On the experimental side, even though 
the d-wave BCS mean field theory captures some features of the superconducting state, the 
recent high resolution data from both ARPES and STM gives a more delicate picture. A new 
feature is the appearance of so called peak/dip/hump structure, which is most clearly seen 
along the [100] direction in the superconducting state for ARPES [Q, and for STM ||. It 
has been suggested j?]] that it stems from the coupling of electrons to the ir resonance observed 
in neutron scattering studies ||. This idea has been further explored both qualitatively^] and 
quantitatively fllPf! , confirming its validity. In addition to the peak/dip/hump structure, there 



are also indications that the quasi-particle peak seems not to be resolution limited[ll]. These 
features indicate the need for a tunneling theory that includes the effect of fluctuations. 
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In this work, we investigate the effects of the spin fluctuations on the tunneling conductance 
spectra along the [100] and [001] (c-axis) directions [12]. By using the Keldysh formulation, we 
first demonstrate that in the high junction resistance limit, the dominant contribution comes 
from the spectral function and the peak/dip/hump structure exists in this limit. This structure 
results from collective spin fluctuations and weakens as the doping 5 increases as reflected in 
some of the recent data, such as in Ref . |0| . Since spin excitations are gapped, they induce 
little qualitative change in the subgap region. Instead, their main effect is to redistribute the 
spectral weight and thus changes the relative strength among currents due to difference tun- 
neling processes. Therefore, to investigate the other limit when the subgap is dominated by 
the Andreev reflection , we define the optimum matching as the condition when the zero bias 
Andreev conductance peak reaches maximum. This corresponds to the Z = case in the BTK 
theory. Under this condition, we show that the Andreev peak can coexist with the dip/hump 
structure, which should be observable in the recent future. Base on our analysis, we will also give 
a possible explanation on a recently observed c-axis STM data |L4|], which shows an unexpected 
step like feature in the negative bias in addition to the peak/dip/hump structure. 

We start by considering a junction consisting of a 2D normal metal on the left (L) hand side 
(— oo < x < 0) and a 2D superconductor (a < x < oo, a is the lattice constant) on the right (R) 
hand side, governed by the Hamiltonian Ht and Hr respectively. The tunneling Hamiltonian 
that connects the surface points at x = and x = a is given by Ht = J2 y t(\yL — Ur\)(c^ l cr + 
c r cl), where the summation is over lattice points along the interface, chosen to be in the y 
direction. We consider the simplest case when the lattice points along the interfaces are equally 
spaced and match the bulk lattice of the metal and the superconductor. The superconductor 
is assumed to have a square lattice with one of the axis parallel to the x direction. The total 
grand Hamiltonian is then given by K = Ht — + Hr — [IrNr + Ht-, where /jll and /ir are 

the chemical potentials and their difference \xl — [J-R is fixed to be the voltage drop eV across 
the junction. 

The tunneling current can be calculated perturbatively by using the Keldysh formalism [fll|. 
This approach was previously applied successfully to study a number of tunneling problems [|i~6|, 



17]. We shall follow Ref. |Tq| and neglect the vertex corrections of Ht- The perturbation series in 



Ht can be then summed exactly. The contribution to the differential conductance G = dl/dV 
can be classified into four terms due to different tunneling processes jjlfj]: G\ is due to particle 
to particle tunneling; G2 is due to particle to particle tunneling with pair creation/annihilation 
as the intermediate state; G3 is due to particle to hole tunneling; and Ga is the Andreev 
conductance. Because Ht is a tight binding model, all the Green's functions have to be replaced 
by the surface Green's functions that connect different points on the junction. Thus there is an 
extra integration over k y and associated with each surface Green's function, there is a t(k y ) factor. 
The function t{k y ) characterizes the spread of the electron wave function along y direction when 
hopping across the junction and can be generally expanded in a cos Fourier series. Note that 
the bare surface Green's function, without being renormalized by Ht, is a 2x2 matrix go(ui, k y ) 
in Nambu's notation[16], and its relation to the bare bulk Green's function Gq(lu, k x , k y ) is given 
by go(uj,k y ) = ^ Jq dk x sin 2 (fe x )Go(o;, k x , k y ), where the momentum k is in unit 1/a. 

The relations of conductance G a to the Green's functions can be best demonstrated in the 
limit when the metal is approximated by its bandwidth ti with a constant density of state, 
i.e., gL,(w,ky) = —i/tz X I, where I is a unit matrix. This avoids complications due to the 
band structure from the metal side. In this case, when t(k y ) = t is a constant, the only 
dimensionless parameter is A = t 2 /(tTtR), where tR is the hopping scale of the superconducting 
side. The junction conductance is then of order A e 2 jfi. For small A, G\(V) is O(A) and is simply 
proportional to the single particle density of state [^qj. Similarly, G^fV) is of order 0(A 2 ) and 

probes / dk y lm(gQ R12 ) , Gs(V) is of the order 0(A 3 ) and probes / dkyPR^i Is'o-R 12 1 2 • Here 1 



and 2 are indices for Nambu notations , Green's functions with the index r are retarded and 

p is the spectral function. Since G3 is subdominant to G2CV) , G2CV) + Gs(V) is negative. 

r -1 2 

Finally, the Andreev conductance is 0(A ) and probes / dk y {g$ R 12 ) • Overall speaking, for 
small A, the total conductance is dominated by G±, corresponding to the large Z limit of the 
BTK theory. In the other limit when Ga dominates in the sub-gap region, the situation is more 
subtle. For s-wave BCS superconductors, the analytic mapping from Z to A was obtained in 
Ref [jig] . The optimum matching (Z = 0) does not occur at large A because the mapping is not 
monotonic. In general, such analytic mapping does not exist , we shall resort to numerics to 
find the optimum matching condition. 

We first consider a simple d-wave BCS superconductor described by 

H R=lZ e k C ka°ka ~ A k( c k^ C -kl + h.C.), (1) 
ka k 

where the dispersion = — 2tji[cos(k x ) + cos(A; y )], and the gap A(k) = Aji[cos(k x ) — cos{k y )]. 
The metal side has the same Hamiltonian with Al = 0. Under the optimum matching condition, 



one obtains a single peak (the Andreev peak) in the total conductance near V = 0[17, 12 1. The 
Andreev peak obtained here is a result of subtle balance between G\ + Ga and G2 + G$. In 
fact, the effect of G2 + G3 is to bring down the quasi-particle peaks in G\(V) so that a single 
peak is manifested. As we shall see, such simple realization of the Andreev peak does not always 
happen in real high-Tc systems due to spin fluctuations. 

To include the spin fluctuations, we shall work with the 2D t — t' — J model. In the slave- 
boson method, the physical electron operators Cj CT are expressed by slave bosons bi carrying the 
charge and fermions representing the spin; Cj CT = bf The mean-field d-wave SC state 

is characterized by the order parameters A = (fi^fji - filfji), Xo = E<x (fi^fja) and the 

condensate of bosons bi — ► (bi) = yS. Eq.(|l]) is then the Hamiltonian for the spinons, /j, with 
dispersion = —2(5tn + J'xo)[cos(k x ) + cos(k y )] — 4:5tR cos(k x ) cos(k y ) — /xr and Ar = 2J'Aq, 
where J' = 3J/8. We shall adopt the following numerical values £r = 2J, tjl = — 0A5tR, and 
J = 0.13eV|lO|, p~9f] . The mean-field parameters xo> an d the chemical potential [1r for different 



doping 5 are obtained from a self-consistent calculation [1C]. Next we include the spin fluctuations 



by perturbing around the mean field Hamiltonian i.e., we write Kr = H R + H', and treat 
H' as a perturbation. In order to account for the tt resonance |8| as well as many other effects of 
spin fluctuations, we calculate the spin susceptibility in a modified random-phase approximation 
(RPA) as defined in Ref. [10, The usual RPA sums over selected sets of graphs for the spin 



susceptibility x as shown in Fig.l(a) and gives rise to x(q, u) = xo(<l) w )/[l + a ^( c l)xo(q] w)] 
with a = 1. Here, J(q) = J(cosq x + cosq y ), xoi^l,^) is the unperturbed spin susceptibility 
due to the spinon bubbles and the tt resonance emerges as the pole of the denominator. In the 
current approach, a is not one and is considered as a phenomenological parameter whose value 
is chosen such that the AF instability occurs right at the experimental observed value 5 = 0.02. 



For the material parameters we adopt, a is 0.34[10|. 



The inelastic scattering of electrons off the spin fluctuations is taken into account by incor- 
porating x i n t° the self energy of the spinons in the lowest-order approximation. In the SC 
state, there are two different self-energies E s and T, w as shown in Fig. 1(b) and (c). The Green's 
function for spinons is calculated by Gf(\s.,u) = [Gj (k, u) + (A/- + Y, w ) 2 GJq(— k, — cj)] _1 with 
G^o(k, uj) = [iuj — ek — £<s(k, a;)] -1 . Since bosons condense, the physical electron Green's function 
can be simply obtained by G(k, u) = <5G/(q, uj), i.e., only the dynamics of spins is considered. 
Following previous prescriptions, one then obtains the surface Green's functions and thus the 
various conductance. The truncation to the lowest order cannot really be justified rigorously 
so far. Its merit rests mainly upon its simplicity and its usefulness in previous applications to 



problems related to spin fluctuations! 10 1. These studies indicate that it has captured the main 



features of ARPES data along [100] direction and for other directions it also reproduces the 



observed cos(6#) deviation from the pure d-wave[10|. Here we shall examine its validity against 
tunneling data. Note also that we had neglected the spatial dependence of the pair potential, 
which is generally considered not important in [100] and [001] directions. 

We first analyze small A limit. Fig. 2 shows the total conductance with RPA correction for 
various dopings. The positions of the peak and hump are seen to scale weakly with doping. When 
doping increases, the height of peak increases with doping, in consistent with experiments [p"3|] , 
at the same time, the width of the peak increases and tends into the hump region so that the 
hump is smeared out in slightly overdoped region. Another feature which can also be observed 
in the data is that the dip/hump feature at positive bias is always weaker. The precise reason 



behind them can be traced back to the underlying structure of e^. In fact, detailed analysis|12] 
shows that the band edge extends to higher positive bias so that the dip/hump is smeared, while 
the band edge for negative bias essentially stays at small bias, leaving the dip/hump unsmeared. 

We now numerically identify the optimum matching condition so that the zero bias Andreev 
conductance peak can be manifested best. For each k y , we compute the optimal value t op t 
such that the Andreev conductance at V = reaches maximum. The resulting t op t(k y ) can be 
approximated by 

t op t(k y ) = a + a\ cos(k y ) + a 2 cos(2/c y ) (2) 

This implies that including next nearest neighbor hopping along the junction is necessary. How- 
ever, the forward hopping ao still dominates (for instance, when 5 = 0.12, we obtain ao = 2.41, 
til = —0.44, and a 2 = 0.34). In Fig. 3, we show the optimal manifestation of the Andreev peak 
for different doping. The metal side is modeled by a simple tight-binding model on the square 
lattice. We see that the dip/hump structure coexists with the Andreev peak. Fig. 4 shows 
a similar plot but now the density of state of the metal side is a constant. In this case, the 
Andreev peak never out wins the quasi-particle peaks resulted from G\{V) so that a plateau 
is observed. In both cases, the trend of the dip/hump structure with doping is consistent with 
what is found in Fig. 2. All these qualitative features should be experimentally verified in the 
future as a test of the mechanism of the spin fluctuations. 

To further test this particular RPA approach, we compute the c-axis tunneling spectrum. 
Fig. 5 shows our numerical results, in comparison to the recent STM curve by Pan et. al. [ [l4[| . 
It is quite encouraging that two curves are very similar in the shape. In particular, the step 
around 45 mV is reproduced in the RPA approach at slightly larger bias. This step results from 
the band edge, which, as we mentioned, essentially stays at small bias as one changes doping. 

To summarize, we have analyzed the effects of spin fluctuations on the SN junction using 
the Keldysh formulation and a modified random phase approximation. The peak/dip/hump 
structure is reproduced and we show that it disappears gradually as one goes to slightly over- 
doped region. Using the same formulation, we predict that the dip/hump structure can coexist 
with the zero bias Andreev peak in optimal matching conditions. Our analysis on the c-axis 
tunneling shows good qualitative agreement between this approach and the recently observed 
STM tunneling curve. 
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1 FIGURE CAPTIONS 



Fig.l Feynman diagrams for (a) the spin susceptibility and (b) and (c) the lowest order contri- 
butions to the self-energy from spin fluctuations. 

Fig. 2 The total conductance with RPA correction in the tunneling limit. Here the metal is 
modelled by a constant density of state with A = 0.05. 

Fig. 3 The optimal manifestation of the Andreev peak with a square lattice (tz = 1-0) for 
the metal side. 

Fig. 4 The optimal manifestation of the Andreev peak with constant density of state (tz = 
1.0) for the metal side. 

Fig. 5 The density of state for c-axis tunneling for 5 = 0.12. Inset: The STM tunneling curve 
observed by Pan et al| 14 ] . 
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